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Free  Internal  gravitational  waves  arise  as  a reaction  of  a stably  stratified 
heterogeneous  liquid  to  an  energetic  external  action.  They  carry  energy  from  the 
region  of  an  initial  disturbance  to  the  ambient  space,  tending  to  return  the  liquid 
to  the  initial  state  of  equilibrium. 

In  the  majority  of  theoretical  studies  of  free  internal  gravitational  waves  in 
the  ocean,  the  basic  assumption  is  that  they  are  small  harmonic  oscillations  of  an 
ideal  incompressible  liquid  about  a certain  state  of  equilibrium.  The  state  of  jl 

rest  is  usually  taken  as  this  state  of  equilibrium.  In  this  case,  the  system  of 
fluid  mechanics  equations  describing  the  unperturbed  state  of  equilibrium  is  reduced 
to  a single  equation  of  statics.  Such  a choice  of  the  principal  state  is  equivalent 
to  neglecting  the  convection  terms  in  linearized  equations  for  perturbations,  in 
comparison  with  local  changes.  V.  Krauss^notes  that  this  operation  is  inadmissible 
in  the  study  of  tidal  internal  waves.  A similar  situation  arises  in  the  study  of 
other  types  of  internal  gravitational  waves,  whose  phase  velocities  are  comparable 
to  the  velocities  of  sea  currents  existing  in  the  initial  state  of  equilibrium. 


Since  the  spectrum  of  phase  velocities  of  internal  gravitational  waves  in  the 
ocean  is  fairly  wide,  and  a certain  average  transport  of  waters  is  the  rule  rather 
than  the  exception  for  all  regions  of  the  World  Ocean,  it  appears  that  sea  currents 
(the  main  flow)  must  be  considered  in  order  to  explain  many  aspects  of  the  dynamics 
of  internal  gravitational  waves.  Sea  currents  may  cause  significant  changes  in  the 
frequencies  of  internal  waves,  inbluding  the  dynamic  instability  of  oscillations. 

The  article  presents  certain  considerations  concerning  the  dependence  of  the 
parameters  of  free  internal  gravitational  waves  on  the  characteristics  of  steady 
sea  currents,  derives  the  necessary  condition  for  the  existence  of  internal  gravi- 
tational waves,  and  explains  the  phenomenon,  frequently  observed  during  actual  ob- 
servations, of  an  abrupt  change  of  phase  of  internal  oscillations  with  depth.  Also  /ll.'  \ 
discussed  is  the  question  of  errors  in  the  determination  of  the  true  periods  of 
internal  waves  based  on  observations  of  hydrological  elements  used  as  indicators 
of  internal  oscillations.  These  errors  are  due  to  the  existence  of  the  main  flow 
and  are  generated  by  the  Doppler  effect.  i 

1.  In  the  case  of  an  incompressible  heterogeneous  liquid,  the  initial  state 
of  which  is  the  state  of  rest,  the  exact  lower  boundary  of  the  periods  of  free 
internal  gravitational  waves  was  first  Indicated  by  Groen.^6  Subsequently,  on  the 
assumption  of  this  principal  state  of  equilibrium,  estimates  of  frequencies  of 
progressive  harmonic  free  internal  gravitational  waves  were  made  in  Refs.  2,  6,  9, 

.10,  15,  and  17.  j 

Q IQ 

Thus,  for  example,  Monin  and  Obukhov  ’ showed  for  an  incompressible  liquid 
that  the  frequencies  of  free  gravitational  waves  1 te  within  the  interval 

2 2 -2 
r $ a < gya  , 


Numbers  in  the  right  margin  indicate  pagination  in  the  original  text. 
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where  l is  the  Coriolis  parameter; 

a is  the  oscillation  frequency; 
g is  the  gravitational  acceleration; 
a is  the  velocity  of  sound; 

da2 

Y = (k-  l)g  + is  the  thermal  stability  parameter; 

K is  the  expansion  exponent; 
z is  the  vertical  coordinate,  directed  upward. 

Since  in  adiabatic  processes  Y = (where  T = -( is  the  Sverdrup- 

\ Pq  “z  • 

Hesselberg  stability  and  is  the  density  in  the  unperturbed  state) , by  introducing 
the  Vaisala-Brunt  frequency  N = ✓gT  and  passing  to  an  incompressible  liquid  (a  ■+•  °°)  , 

we  find  that  in  the  ocean,  the  frequencies  of  the  waves  in  question  satisfy  the 

2 2 2 
inequality  i 4 0 < N . 

Reference  2 discusses  the  case  of  an  isothermal  atmosphere  in  the  absence  of 
the  Coriolis  force  and  shows  that  all  the  frequencies  of  free  gravitational  oscilla- 
tions are  bounded  by  the  frequency 

aQ  = /(.<  - l)gK~1h~1, 

where  h is  the  height  of  the  atmosphere. 


In  the  general  case,  the  height  of  the  equivalent  homogeneous  atmosphere  is 

2 -1  -1 

related  to  the  sound  velocity  a by  the  relation  h(z)  = a K g , and  for  an  isothermal 
atmosphere  Y = ~ l)g* 


Obviously,  in  an  incompressible  liquid,  this  case  corresponds  to  a constant 
gradient  of  the  density  logarithm,  independently  of  the  height,  and  the  set  of  gra- 
vitational wave  frequencies  has  the  following  frequency  as  the  upper  bound: 


N = /- g In  Pq  = const. 

It  will  be  shown  that  the  characteristics  of  internal  gravitational  waves,  / 
treated  as  small  oscillations,  should  depend  on  the  properties  of  the  initial 
state  of  equilibrium. 

Let  the  sea  be  an  ideal,  incompressible,  stably  stratified,  horizontally  in- 
finite liquid  of  thickness  H.  We  will  assume  that  the  density  of  this  layer  Pq  as 
well  as  the  velocity  ug  and  direction  of  horizontal  mass  transfer  therein  are  func- 
tions of  only  the  vertical  coordinate  z.  We  direct  the  z axis  of  the  rectangular 
coordinate  system  vertically  upward,  and  the  x and  y axes,  horizontally.  We  place 
the  origin  at  a flat  horizontal  seabed.  We  denote  by  8(z)  the  angle  between  the 
positive  direction  of  the  x axis  and  the  direction  of  mass  transfer  at  any  fixed 
level  in  the  layer.  Let  the  perturbation  in  the  liquid  be  caused  by  a progressive 
free  internal  gravitational  wave  traveling  at  angle  a to  the  positive  direction  of 
the  x axis.  We  represent  the  internal  oscillations  as  small  harmonic  perturbations 
of  the  three  velocity  components,  pressure  and  density  in  relation  to  the  principal 
state  of  equilibrium.  Let 


r'  = r(z)exp[i(at  - kx  - sy) ] , 
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where  r'  is  any  of  five  perturbed  hydrodynamic  elements; 
r is  a complex  amplitude  factor; 

k and  s are  horizontal  wavenumber  components  for  which  a * arctan  s/k; 
t is  the  time; 
i is  an  Imaginary  unit. 


We  then  have  the  following  second-order  ordinary  homogeneous  linear  differen- 
tial equation  in  the  amplitude  of  the  vertical  velocity  of  internal  waves  w:* 

Ip { r ^ dv  . 7 du  1)  dv 

{ P—  m7{c  — Ju)»  |_  4z  {c  — u)  ' dz  If  ~dz 

— 77 — =—  {a*  — (c  — it)  \m'  (c  — It) -£-(o0  —)]  — 

/— «i'((-ii)!  l l ' f0  dz  V 0 dz  I) 

^ f J HLfn  dv  ^ i 1 4 u dv  11  ... 

m l tv  dz  V®  dz)  (c  — u)  ' dz  ' -dz  Jjw  — 0. 


where 


u = //0  cos  (a  — |J); 
v = u0  sin  (a-?); 
c - °m~\ 

The  boundary  conditions  for  Eq.  (1)  will  be  written  as 

w(0)  = w(H)  = 0 


(2) 


With  l = 0,  we  obtain  from  expression  (1) 
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0.  (3) 


We  will  hereinafter  assume  that  the  function  u has  continuous  derivatives  up  to 
second  order  inclusive,  and  the  function  Pq  has  a continuous  first  derivative  in- 
side the  layer  [0,  H]. 

The  conclusion  that  there  is  a lower  period  of  internal  waves  in  the  presence 
of  a steady  main  flow  is  obtained  most  simply  from  the  equation  for  the  deviation 
£ of  the  liquid  particles  from  the  equilibrium  position,  related  to  the  amplitude 
of  the  vertical  velocity  w as  follows: 

w(z)  = £(z)im(c  - u).  (4) 


From  expressions  (3)  and  (4),  we  obtain  an  equation  in  £: 


• 

d7; 
dz 7 


r -| — — . ~ 

(C  — //  ) lfz 


+ f—  A-~ w'1;  = 0. 

* [ (c  — u)*  J 


(5) 


Equation  (5)  should  be  examined  for  homogeneous  boundary  conditions  £(0)  = £(H)  = 0 


This  equation  can  be  obtained  by  simple  transformations  of  the  equation  given 
in  the  monograph  of  V.  Krauss.^ 
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In  deriving  inequality  (6) , c was  assumed  to  be  a real  quantity  having  the 
meaning  of  phase  velocity  of  the  internal  wave. 


This  method  was  used  by  M.  Yasui^O  in  a study  of  the  stability  conditions  of 
internal  waves  traveling  in  the  frontal  zone  of  a horizontally  inhomogeneous  ocean. 


To  estimate  the  lower  period  of  free  internal  gravitational  waves,  we  will  use 
Groen's*  method,16  applying  it  to  Eq.  (5).  Since  the  function  £ inside  the  interval 
is  assumed  not  identically  equal  to  zero  and  twice  continuously  differentiable,  for 
the  stated  boundary  conditions,  according  to  the  Rolle  theorem,  it  should  reach  an 
extreme  value  in  this  interval.  In  this  case,  there  is  at  least  one  level  z.  for 
which  J 


Hence,  according  to  Eq.  (5),  the  following  inequality  holds  for  such  levels: 

N2(c  - u)"2  - m2  > 0.  (6) 


Considering  that  the  internal  wave  period  T 
lowing  estimate  from  inequality  (6): 

’>!ri«-4|. 


= 2Tr(mc)  , we  obtain  the  fol- 


Hcnce,  the  following  inequality  is  all  the  more  valid: 
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where  Tq  = 27tNq^  is  the  Vaisala  period,  where  -^o  Therefore,  in  the 


presence  of  the  main  flow,  there  exists  no  accurate  lower  period  of  free  internal 
gravitational  waves.  The  actual  period  of  these  waves  may  be  larger  or  smaller 
than  the  Vaisala  period,  depending  on  the  sign  of  the  projection  of  the  velocity 
of  steady  flow  on  the  direction  of  propagation  of  the  wave,  with  the  period  of 
internal  waves  X ->  +0  in  the  coincidence  layer  (c  = U).  The  meaning  of  the  latter 
statement  is  consistent  with  Krauss'  conclusion6  that  internal  waves  disappear  at 
the  critical  depth,  i.e.,  the  level  corresponding  to  the  coincidence  layer. 


2.  Applying  to  Eq.  (3)  arguments  similar  to  those  given  in  the  derivation  of 
inequality  (6) , we  arrive  at  the  inequality 


, A'1 

rn-  < — 

(c  — //  )* 


Obviously,  inequality  (7)  can  be  fulfilled  only  when 


Hence,  internal  waves  with  a given  phase  velocity  in  the  presence  of  steady 
flows  in  the  liquid  can  exist  only  when  there  exists  at  least  one  layer  for  which 
inequality  (8)  Is  fulfilled.  Let  us  note  that  when  u = const  (in  the  special  case 
TT  = 0) , a known  condition  of  the  existence  of  internal  waves  follows  from  inequality  (8) 


namely:  N >0,  i.e. , If  the  state  of  uniform  motion  (or  the  state  of  rest)  is  the 
principal  state  of  equilibrium,  internal  waves  can  exist  only  when  the  inhomogeneous 
liquid  under  study  includes  layers  in  which  the  density  decreases  with  height. 

3.  We  will  examine  the  question  of  limiting  frequencies  of  free  internal  gra- 
vitational waves  from  the  standpoint  of  the  eigenvalues  of  Eq.  (1).  We  will  assume 
for  simplicity  that  a = 6 = 0 (plane  motion).  Then  Eq.  (1)  will  be  written  in  the 
form 


<p'i> 

_rr 

/»  du^ 

dw 

m* 

dz'- 

l 

(f  -«„)[/>-  m'-  (C  -«„)’]  dz] 

dz 

/*  - »>’  (c  - II')'- 

X (/V2  - (c  - //,)  [m2  (<r - h0)  - (?o  ■$)][ = 0. 

To  obtain  the  dispersion  equation,  we  set 

w = w sin  8z,  where  w = const, 
n n 

Then,  boundary  conditions  (2)  will  be  met  when 

nir 


(9) 


(10) 
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where 


n = 1,  2...  is  the  mode  of  the  internal  wave; 

0 = 0(n)  is  the  vertical  wavenumber. 

Substituting  equality  (10)  into  Eq.  (9),  we  obtain 

0 - sin  rJZ  -f-  / 1'  — n ~rr)  5 cos  1)2  -r 


/»  — m-  (c  — 


-,X 


(c  — «o)  [I7  — m‘‘  \c  — Il0y-'|  dz  f 

X [at  - nf  (c  - uj1  + 1 (c  - H0)  JL  (o0  ^-)]  sin  Hz  = 0 . (It) 

We  will  consider  the  levels  z = z*j  for  which  cos  8z  = 0,  i.e.,  z*^  = Hn-1(j  + 

Since  we  are  interested  only  in  z*.£  (0,  H) , for  each  fixed  n we  have  j = 0,  1,  2..., 
n-1 . J 

For  levels  z*j , Eq.  (11)  will  be  written 

mV  (l  - (0! m-)  0 -1'-  i-  n-  [.V2  -r  P0'  (c  - «,)  X 

**-(*$)]• 


(12) 


Hence 

where  A **  2Trm-^  is  the  length  of  the  internal  wave. 

Ttie  right-hand  side  of  equality  (12)  is  positive  in  accordance  with  inequality  (8). 


If  Uq  H 0,  we  find  with  a high  degree  of  accuracy  from  expression  (12)  that 
the  frequency  of  fairly  long  waves  O l,  i.e.,  approaches  the  frequency  of  iner- 
tial oscillations,  and  the  frequency  of  very  short  waves  0 -*■  N,  i.e.,  approaches 
the  Vaisala-Brunt  frequency. 

In  the  absence  of  currents  in  a horizontally  homogeneous  ocean,  the  frequency 
of  internal  gravitational  waves  is  invariant  relative  to  their  direction  of  propaga- 
tion. It  is  easy  to  see  that  in  the  presence  of  a horizontal  current,  it  is  sub- 
stantially dependent  on  the  direction  of  the  current,  its  velocity,  and  the  direc- 
tion of  propagation  of  the  wave.  Indeed,  in  the  presence  of  the  main  flow,  for  /12( 

fairly  long  waves,  we  find  from  equality  (12) 


1 __ 

1 C ,* 


and  for  very  short  waves 


«’  [,V*  + f;-'  (c  - «,)  £ Tr);  - T,r  • 

<-  2n]  1 - (<-".)  T,  (p*  Tr)]  ”•  (H) 


It  is  significant  that  in  general,  in  the  presence  of  the  main  flow,  there  may 
exist  free  internal  gravitational  waves  with  periods  exceeding  the  period  of  inertial 
oscillations.  As  follows  from  formula  (13),  this  takes  place  when  a very  long  wave 
travels  in  the  direction  opposite  to  that  of  the  main  flow,  this  excess  being  directly 
proportional  to  the  ratio  of  the  velocity  of  the  main  flow  to  the  velocity  of  the 
internal  wave. 


All  of  the  above  considerations  pertain  to  mutually  stable  internal  waves  (c 
being  a real  quantity).  Consideration  of  the  possible  dynamic  instability  of  internal 
oscillations  apparently  will  not  lead  to  any  fundamentally  different  conclusions  con- 
cerning the  effect  of  currents  on  the  frequency  of  internal  oscillations,  but  may  sig- 
nificantly limit  the  actual  frequency  region  of  existence  of  internal  waves. 

Let  us  note  that  the  direct  application  of  Groen's  method  to  Kq.  (9)  will  lead 
to  the  double  inequality 


2*1 1 - 


•_*!! ! i a*. 

77.7  x < ' < / . I 1 r ’ 


which  confirms  the  validity  of  the  approximations  used  in  the  derivations  of  formulas 
(13)  and  (14). 

4.  The  literature  frequently  metions  the  phenomenon  of  oscillation  phase  dis- 
continuity in  actual  observations  of  internal  waves. 3,4,11,13  The  phase  discontinuity 
of  internal  oscillations  is  easy  to  explain  if  it  is  observed  in  layers  of  density 
discontinuity,  since  the  conclusions  of  interface  wave  theory  can  be  applied  to  this 

r.'lQO  . 7 
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In  a continuously  stratified  liquid,  the  presence  of  a phase  discontinuity  of 
internal  oscillations  can,  generally  speaking,  be  formally  explained  on  the  basis 
of  a determination  of  the  normal  oscillations  (modes)  of  the  system.  In  this  case, 
it  can  easily  be  shown  that  the  depth  at  which  the  phase  discontinuity  is  located 
in  disturbed  motion  does  not  necessarily  coincide  with  the  level  of  the  density 
discontinuity  layer.  Moreover,  the  level  at  which  the  phase  discontinuity  takes 

place  can  be  accurately  determined  in  this  case  as  a function  of  the  mode  of  the 
internal  wave. 
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It  is  our  view  that  the  phenomenon  of  phase  discontinuity  may  also  be  explained 
by  using  the  theory  of  hydrodynamic  stability  of  the  motion  of  a heterogeneous 
liquid. 

Tollmin®  showed  theoretically  in  the  case  of  the  motion  of  a viscous  homogeneous 
liquid  that  a phase  discontinuity  in  oscillations  determining  the  perturbing  motion 
may  arise  during  the  passage  through  the  critical  layer.  This  conclusion  was  con- 
firmed by  the  experiments  of  Schubauer  and  Scramsted. 

For  a heterogeneous  ideal  liquid,  a shift  in  the  horizontal  velocity  of  the 
main  flow  frequently  causes  a phase  discontinuity  during  the  passage  through  the 
critical  layer.  This  assertion  results  from  an  analysis  of  the  work  of  Miles. 18,19 

It  may  be  postulated  that  the  phase  discontinuity  for  a given  mode  of  an  internal 
wave  will  occur  at  the  level  where  the  phase  velocity  of  this  mode  is  equal  to  the 
projection  of  the  velocity  of  the  main  flow  on  the  direction  of  propagation  of  the 
wave.  This  assumption  admits  in  principle  an  experimental  check  involving  special 
observations  that  may  be  carried  out  by  means  of  standard  oceanographic  instruments. 

Let  us  note  that  a similar  phase  discontinuity  phenomenon  during  the  passage 
through  the  coincidence  layer  takes  place  in  the  process  of  generation  of  wind- 
generated waves. 14 

5.  In  oceanographic  observations,  it  is  sometimes  useful  to  differentiate  the 
transformation  of  the  wave  process  by  currents,  i.e.,  to  consider  the  Doppler  effect. 
In  a study  of  internal  tidal  waves,  such  an  estimate  was  made  by  Bukhteyev  on  the 
assumption  that  the  waves  travel  along  the  phase  boundary  of  two  homogeneous  liquids 
with  a uniform  principal  motion  in  the  entire  system.  These  results  can  be  easily 
extended  to  the  case  of  internal  gravitational  waves  of  any  period,  traveling  in 
an  arbitrary  direction  in  a continuously  stratified  inhomogeneous  liquid,  in  which 
the  horizontal  velocity  vector  of  the  main  flow  is  an  arbitrary  function  of  the 
vertical  coordinate. 

Then  the  distortion  AT  of  the  period  of  internal  waves  (difference  between  the 
observed  T^  and  true  T2  periods)  may  be  obtained  from  the  equality 


M c -11  1 / 


(15) 


Since  in  this  case  u = u(z),  the  magnitude  of  the  distortion  of  the  internal  wave 
period  substantially  depends  on  the  relation  between  the  phase  velocity  of  the  wave 
and  the  value  of  the  projection  of  the  velocity  of  the  main  flow  on  the  direction  of 
propagation  of  the  wave.  Hence,  in  the  general  case,  the  values  of  AT  will  be  dif- 
ferent at  different  observation  levels. 


/I  2. 
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The  Doppler  effect  is  usually  considered  when  processing  observations  made 
from  a moving  ship.  The  corresponding  procedure  was  given  by  Krauss-*  and  Sabinin.12 
However,  in  processing  observations  of  internal  waves,  made  from  an  anchored  ship 
or  by  means  of  instruments  mounted  on  autonomous  buoy  stations,  the  Doppler  effect 
is  usually  neglected.  It  should  be  kept  in  mind,  however,  that  the  internal  wave 
spectra  obtained  may  be  markedly  distorted  because  the  velocity  of  the  sea  currents 
at  the  level  of  the  observations  is  not  zero. 


If  the  actual  observations  of  internal  waves  are  organized  so  that  one  can  cal- 
culate the  phase  velocity  of  a wave,  direction  of  its  propagation,  as  well  as  the 
speed  and  direction  of  nonperiodic  currents,  then  according  to  equality  (15),  the 
true  wave  period  can  be  determined  from  the  formula 


i *•  i 


We  will  note  in  conclusion  that  the  "observational"  Doppler  effect  and  fre- 
quency change  effect  (discussed  in  sections  1 and  3),  physically  caused  by  sea 
currents,  are  not,  generally  speaking,  identical.  In  particular,  it  appears  that 
resonance  phenomena  in  the  ocean  occur  at  frequencies  modified  by  the  presence  of 
the  main  flow. 
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